Gravitational collapse in f(R) theories 
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We study the gravitational collapse in modified gravitational theories. In particular, we analyze a 
general f(R) model with uniformly collapsing cloud of self-gravitating dust particles. This analysis 
shares analogies with the formation of large-scale structures in the early Universe and with the 
formation of stars in a molecular cloud experiencing gravitational collapse. In the same way, this 
investigation can be used as a first approximation to the modification that stellar objects can suffer 
in these modified theories of gravity. We study concrete examples, and find that the analysis of 
gravitational collapse is an important tool to constrain models that present late-time cosmological 
acceleration. 

PACS numbers: 98.80.-k, 04.50.+h 
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INTRODUCTION 

In the general study of astrophysical weak gravita- 
tional fields, relativistic effects tend to be ignored. How- 
ever, there are clear examples of stellar objects in which 
these effects may have important consequences, such as 
neutron stars, white dwarfs, supermassive stars or black 
holes. Indeed, it becomes necessary to consider observa- 
tionally consistent gravitational theories to study these 
objects. General Relativity (GR) has been the most 
widely used theory but other gravitational theories may 
be studied for a better understanding of the features and 
properties of such objects and to compare their predic- 
tions with experimental results. 

The gravitational collapse for a spherically symmet- 
ric stellar object has been extensively studied in the GR 
framework (see [l[ and references therein). By assum- 
ing the metric of the space-time to be spherically sym- 
metric and that the collapsing fluid is pressureless, the 
found metric interior to the object turns to be Robertson- 
Walker type with a parameter playing the role of spatial 
curvature and proportional to initial density. The time 
lapse and the size of the object are given by a cycloid 
parametric equation with an angle parameter tp. Further 
results are that the time when the object gets zero size 
is finite and inversely proportional to the square root of 
the initial density. Finally, the redshift seen by an exter- 
nal observer is nevertheless infinite when time approaches 
the collapse time. 

In spite of the fact that GR has been one of the most 
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successful theories of the twentieth century, it does not 
give a satisfactory explanation to some of the latest cos- 
mological and astrophysical observations with usual mat- 
ter sources. In the first place, a dark energy contribution 
needs to be considered to provide cosmological accelera- 
tion whereas the baryonic matter content has to be sup- 
plemented by a dark matter (DM) component to give 
a satisfactory description of large scale structures, rota- 
tional speeds of galaxies, orbital velocities of galaxies in 
clusters, gravitational lensing of background objects by 
galaxy clusters, such as the Bullet Cluster, and the tem- 
perature distribution of hot gas in galaxies and clusters 
of galaxies. All these evidences have revealed the inter- 
est to study alternative cosmological theories. This extra 
DM component is required to account for about 20% of 
the energy content of our Universe. Although there are 
many possible origins for this component Q , DM is usu- 
ally assumed to be in the form of thermal relics that 
naturally freeze-out with the right abundance in many 
extensions of the standard model of particles Fu- 
ture experiments will be able to discriminate among the 
large number of candidates and models, such as direct 
and indirect detection designed explicitly for their search 
[U, or even at high energy colliders, where they could be 
produced 

A larger number of possibilities can be found in the 
literature to generating the present accelerated expan- 
sion of the Universe [6j. One of these methods consists 
of modifying Einstein's gravity itself 0, HJ without in- 
voking the presence of any exotic dark energy among the 
cosmological components. In this context, functions of 
the scalar curvature when included in the gravitational 
action give rise to the so-called f(R) theories of modi- 
fied gravity [l]. They amount to modifying the l.h.s. of 
the corresponding equations of motion and provide a geo- 
metrical origin to the accelerated cosmological expansion. 
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Although such theories are able to describe the acceler- 
ated expansion on cosmological scales correctly [l(J , they 
typically give rise to strong effects on smaller scales. In 
any case, viable models can be constructed to be com- 
patible with _local gravity tests and other cosmological 
constraints 



in. 



The study of alternative theories of gravitation requires 
establishing methods able to confirm or discard their va- 
lidity by studying the cosmological evolution, the growing 
of cosmological perturbation and, at astrophysical level, 
the existence of objects predicted by GR such as black 
holes 12] or dust clouds forming compact structures. It 
is well-known that f(R) gravity theories may mimic any 
cosmological evolution by choosing adequate f(R) mod- 
els, in particular that of ACDM [lCf. This is the so-called 
degeneracy problem that some modified gravity theories 
present: accordingly, the exclusive use of observations 



such as high-redshift Hubble diagrams from SNIa 13 1, 
baryon acoustic oscillations or CMB shift factor 15 1 , 
based on different distance measurements which are sen- 
sitive only to the expansion history, cannot settle the 
question of the nature of dark energy [l6| since identi- 
cal results may be explained by several theories. Nev- 
ertheless, it has been proved that f(R) theories - even 
mimicking the standard cosmological expansion - provide 
different results from ACDM if the scalar cosmological 
perturbations are studied [17[. Consequently, the power 
spectra would be distinguishable from that predicted by 
ACDM 0. 

It is therefore of particular interest, to establish the 
predictions of f(R) theories concerning the gravitational 
collapse, and in particular collapse times, for different 
astrophysical objects. Collapse properties may be either 
exclusive for Einstein's gravity or intrinsic to any covari- 
ant gravitation theory. On the other hand, obtained re- 
sults may be shed some light about the models viability 
and be useful to discard models in disagreement with ex- 
pected physical results. 

In [lj| the authors studied gravitational collapse of a 
spherically symmetric perfect fluid in f(R) gravity. By 
proceeding in a similar way to [![, the object mass was 
deduced from the junction conditions for interior and ex- 
terior metric tensors. Finally, they concluded that /(i?o) 
(constant scalar curvature term) slows down the collapse 
of matter and plays the role of a cosmological constant. 
Authors in (20| paid attention to the curvature singular- 
ity appearing in the star collapse process in f(R) the- 
ories. This singularity was claimed to be generated in 
viable f(R) gravity and can be avoided by adding a R a 
term. They also studied exponential gravity and the time 
scale of the singularity appearance in that model. It was 
shown that in case of star collapse, this time scale is much 
shorter than the a ge o f the universe. Analogous studies 
were carried on by |21| claiming that in this class of the- 
ories, explosive phenomena in a finite time may appear 
in systems with time dependent increasing mass density. 



Reference [22j includes a complete study of neutron 
stars in f(R) theories. The most relevant result in this 
investigation suggests that f(R) theory allows stars in 
equilibrium with arbitrary baryon number, no matter 
how large they are. Very recently authors in [23| studied 
collapse of charged black holes by using the double-null 
formalism. 

Charged black holes in f(R) gravity can have a new 
type of singularity due to higher curvature corrections, 
the so-called f(R)-induced singularity, although it is 
highly model-dependent. 

The present work has been arranged as follows: in sec- 
tion II, f(R) modified gravity theories will be introduced. 
Gravitational collapse in f(R) theories will be presented 
in section III. After performing some calculations, the 
evolution equation for the object scale factor will be ob- 
tained. This equation will be used throughout the fol- 
lowing sections. Section IV is then dedicated to achieve 
solutions for the modified equations in three qualitatively 
different f(R) models, which try to illustrate the broad 
phenomenology of the subject. This is therefore the aim 
of this section: to study gravitational collapse by calcu- 
lating the evolution of the object scale factor in particu- 
lar f(R) models. Finally, the conclusions based upon the 
presented results will be analyzed in detail in section V. 



II. f(R) THEORIES OF GRAVITY 

With the aim of proposing and alternative theory to 
GR, a possible modification consists of adding a func- 
tion of the scalar curvature, f(R), to the Einstein- Hilbert 
(EH) Lagrangian. Therefore the gravitational action be- 
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d?xy/\g\(R + f(R)) 



(1) 



By performing variations with respect to the metric, the 
modified Einstein equations turn out to be 

(I + fn)R fW - ^{R + f(R))g^ + V^ u f R = -8nGT^, 

(2) 

where T^ v is the energy-momentum tensor of the matter 
content, Jr = df(R)/dR and V^ v = V M V„ — g^D with 
□ = V Q V" and V is the usual covariant derivative. 

These equations may be written a la Einstein by iso- 
lating on the l.h.s. the Einstein tensor and the f(R) 



1 In the present work we employ the natural units system in which 
h = c = 1. Note also that our definition for the Riemann tensor 
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contribution on the r.h.s. as follows 

Rau — t^RQij-v = y\ — I — FT — 87rGT„y 
1 I 1 + JR) 



(3) 



We can also find the expression for the scalar curvature 
by contracting ([2]) with g^ v which gives: 



(1 - f R )R + 2f(R) + 3Df R = SttGT. 



(4) 



Note that, unlike GR where R and T are related alge- 
braically, for a general f(R) those two quantities are 
dynamically related. In the homogeneous and isotropic 
case, the scalar curvature in f(R) theories becomes 



R 



87rGT-2f(R)-3f R 



(5) 



where dot means the derivative with respect to cosmic 
time. 



III. GRAVITATIONAL COLLAPSE IN f(R) 

In the case of our investigation, we introduce the spher- 
ically symmetric metric 

ds 2 = dt 2 - U{r, t)dr 2 - V(r, t)(d6 2 + sm 2 8dcj) 2 ) (6) 

If the collapsing object is approximated to be pressureless 
p ~ 0, the components of the energy- momentum tensor 
can be expressed as follows 

= pu^Uy ; T t t = p; T\ = if i = r,9, (j). (7) 

We may further simplify the collapse model by consider- 
ing p independent from the position. Therefore, we can 
search -as is actually the usual approach in the GR case- 
a separable solution for this metric as follows 



U(r,t) = A 2 (t)h(r), V(r,t) = A 2 (t)r 2 



(8) 



where a previous reparametrization of the radial coor- 
dinate is required. When f(R) modified tensorial equa- 
tions are studied in the homogeneous and isotropic case 
- in which f(R) does not depend on the position-, the 
trace component provides 
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From ((HJ), we deduce that A\ and A 2 are proportional, in 
other words, Aj(t) = C(r)A 2 (t). So, if we choose Ai(t) — 
A 2 (t) = A(t), the dependence in the radial coordinate is 
reabsorbed by h(r). Hence: 



U(r,t) = A 2 (t)h(r), V(r,t) = A 2 (t)r 2 



(10) 



Components tt, rr and 88 for the modified tensorial equa- 
tions may be written respectively in terms of the func- 
tions A{t) and h{r) as follows 



o A - 1 

A (1 + f R ) 



-8TrGp + 3^f R 



~(R + f(R)) 



(11) 
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(12) 
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(R + f(R)) 



(13) 



Let us point out two important aspects of equations (|12[) 
and (fT3| : firstly, terms on the r.h.s of both equations are 
equal. Secondly, the term on the l.h.s. exclusively de- 
pends on r, whereas the term on the r.h.s. only depends 
on t in both equations, so that they must be constants 2 . 
Therefore we may equal l.h.s. of both equations to pro- 
vide 
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(14) 



where we have equaled both equations (multiplied by a 
factor A 2 ) to a constant —2k. The resulting solution is 
h(r) = (1 — kr 2 )^ 1 . Once we have calculated h(r), the 
resulting metric can be expressed as follows: 



ds 2 



dt 2 -A 2 {t) 



dr 2 



1 - kr 2 



\d8 2 



. (15) 

which is formally the same as the one obtained in the 
GR case [l| . Expression (fT4|) for k may be substituted in 
either expression (fT2|) or (fl"3)) yielding: 



A 
A 



- 2 



2k _ 1 
A* ~ (1 + Jr~) 



-fR ~ ^~^fR 



f{R)) 



(16) 



Taking into account p(t) = p(t = 0)/A(t) 3 (given by the 
energy motion equation for dust matter) and the results 



2 This fact is also satisfied in the GR case and allows the simplifi- 
cation of the calculus. 
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in (fT6)l . equation (fTTj) becomes: 
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-^(O)^- 1 + -A 2 f R 



(R + f(R)) . (17) 



Furthermore, provided that the fluid is assumed to be 
at rest for t = 0, initial conditions A(t = 0) = and 
A(t = 0) = 1 hold. This last condition means that the 
scale factor of the object at initial time is normalized 



to unity. In order to simplify the notation we define 
R(t = 0) = i?o and p(t — 0) = p . Therefore, evaluation 
of (fl7|) at t — allows to recast k as follows 
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f(Ro)) 

(18) 



Once k has been expressed in terms of different quantities 
initial values, equations (J5J) and (|18[) may be inserted in 
(|T7| to provide 
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6(1 - p R (Ro)) 
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(19) 



The previous expression will be solved perturbatively to 
first order in perturbations for different f(R) models. Let 
us remind at this stage that the zeroth order solution of 
GR is given by the parametric equations of a cycloid [l[ : 



ip + sin -0 
2Vk ' 



A G = -(l + cos^). 



(20) 



Expression (|20j) clearly implies that a sphere with initial 
density po and negligible pressure will collapse from rest 
to a state of infinite proper energy density in a finite time 
that we will denote Tq . This time is obtained for the first 
value of ip such as Aq = 0, i.e. for ip = ir. It means 



T G = 



2Vk 



7T / 



2Vk 2 \8irGpo 



1/2 



(21) 



In order to study the modification to the gravitational 
collapse in f(R) theories, we will expand A around Aq 
and f(R) around the scalar curvature in GR (R = Rg)'- 

A = Aq + g(ip) , (22) 



f(R)~f(RG) + f'(RG)(R-R G ). (23) 

The presence of a function f(R) in the gravitational La- 
grangian will represent a correction of first order with 
respect to the usual EH Lagrangian. Hence, g(ip) as de- 
fined in |22|) will be also first order at least. By substitut- 
ing the series expansions (|22|) and (|23|) in expression (fT9|) 
until first order in e, we find that equation (|19|) becomes 



I 

where we have cancelled out the GR exact solution and IV. f(R) THEORY RESULTS 

only kept first order perturbed terms. Equation (|24l) will 
provide g(ip) evolution for different f(R) models to be 
considered in the next section. 



In this section we shall consider three illustrative f(R) 
models and study the gravitational collapse process for 
collapsing dust. The models under consideration are 
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Model 1: f(R) = eR 2 

This function has been proposed both as a viable in- 
flation candidate [Ht and as a dark matter model [jjj ]. 
In this last reference, the e parameter definition reads 



1 



6toq 



2 " 



(25) 



and the minimum value allowed for mo is computed as 
mo = 2.7 x 10~ 12 GeV at 95 % confidence level, i.e. 
e < 2.3 x 10 22 GeV~ 2 . On the other hand, e > is 
needed to ensure the stability of the model, since in the 
opposite case, a tachyon is present in the theory. These 
constraints are in agreement with 26]. 



After some algebra, for this model, equation (l24l) can 
be written as follows: 



g'(i/>) + g(i>) csc(^) + ^ke (sin(^) + 2 tan ( | 



4<au' f | ) <><•(,•) j = 0. 



(26) 



The homogeneous equation associated to (|26|) presents 
the solution ghomiib) oc cot(V>/2), which diverges at 
(/> = 0. Therefore, its contribution will be ignored in 
the upcoming analysis. The analytical full solution of 
(|2"6")) becomes 



Cl cot ( — 



4^-ot(t 

128 V 2 



64 fib 
16 + sin(^)) + — - tan I ^ 



sec 2 f I) (sec 2 (| 



(27) 



This analytical solution given by (127)) can be compared 
with the GR one by plotting them together as shown in 
Figure [TJ As we see in this Figure, in the first stage of 
the collapse, the correction is negative, what implies that 
we have a larger contraction. On the contrary, very close 
to ip = 7r, where the solution can be approximated as 



72ke 



5(^-tt) 4 



(28) 



the sign of the modification changes and the total collapse 
is avoided. Exactly at this moment, the perturbation 
leaves the linear regime and a more complete analysis 
is required. It is interesting to estimate when the linear 
approach fails and an important modification is expected. 

By using the collapse time parametrization (120[) in the 
ib — ^ 7r limit, one gets 



t 



ib + sin(V>) 7T + 1 /Q(Aibf 



2y/k 
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or written in terms of the relative variation: 



At 

tGR 



67T 



(29) 



(30) 



We are interested in estimating the region of the pa- 
rameter space of the model, where the modified col- 
lapse, and the result for %bc {^P value for the collapse) 
is significantly different from the one predicted by GR 
("0c,Gi?. = k). With this purpose, we can estimate the 
values for which Aq is of the same order of its correction. 
As one can see in Figure [TJ this deviation is more impor- 
tant close to the final stage of the collapse, for ip ~ ir. In 



this region, Aq can be approximated by: 

l , , (ib — 

A G = -(l + cosV)^ 4 



(31) 



We can use these approximations in the limit ib — > tt to 
determine the intersection between the particular solu- 
tion of (|26|) and the GR solution: Aq- This calculation 
will help establishing the validity regime of the pertur- 
bative approach. Therefore, by imposing |g('0)| = |^4g|j 
with g(ib) given by equation (1281) . we obtain: 



ib 



/288|fce| 
V 5 



1/6 



(32) 



At this point, it is necessary to clarify the physical 
value for k in order to discuss if the departure from lin- 
earity is important, k is the initial condition given by 
equation (fl"8|) that depends on the matter density, the 
initial curvature and the particular f(R) model. In our 
analysis we are interested in studying the modification 
to the gravitational collapse in GR and for this reason 
we will assume the same value of k than as given in GR. 
This implies that the entire modification has a dynam- 
ical origin and it does not come from a change in the 
initial conditions. Therefore, we will assume that k only 
depends on the matter density: 



k 



IttG 
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For the most physically interesting values of k and e, 
for which we have studied gravitational collapse of a dust 
matter cloud, the value of ib is quite close to ib = tt and 
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Figure 1: Comparison between the solution given by ||27|I setting 
ci = and the GR case for the model 1: f(R) = eR 2 . The plotted 
k value is fixed by Eq. H33I I. where the density is psF — 1-5 X 10~ 38 
GeV 4 ~ 3.5 X 10 — 18 kg/m 3 , i.e., the matter density in the early 
Universe at redshift z ~ 1100 marking the decoupling of matter 
and radiation and the beginning of structure formation (SF). The 
modification is extraordinarily small and has been increased 52 
orders of magnitude to make it observable: g(ip) = 10 52 g(t/>). 



therefore the asymptotic approach to obtain (|32|) is fully 
justified. The results are summarized in Figure [2j were 
the non- linear regime is shown for different values of e and 
initial densities. For example, it is interesting to check 
the behavior for the matter density in the early Universe 
at redshift z ~ 1100, which marks the decoupling of mat- 
ter and radiation and the beginning of structure forma- 
tion (SF): psF ^ 1-5 x 10- 38 GeV 4 ~ 3.5 x 10- 18 kg/m 3 . 
In this particular case, the calculated root presents a 
slight difference with respect to the solution for GR. In 
particular, ipe differs from ipe, GR = tt in the ninth sig- 
nificant figure. 

Although this modification is not detectable for this 
first model as we deduce from the last considerations, 
it is interesting to stress that the relative modification 
is higher for denser media since the correction increases 
with po as At /ten oc y/po- This behavior is significantly 
different with respect to other models as we will see in the 



following sections. In Figure [2] we have represented some 
relevant density values as well as the non linear regime 
region. The loss of the linear regime takes place at high 
densities since the correction is directly proportional to 
P- 
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Figure 2: Validity of the perturbative regime for model 1, showing 
different relevant regions: In blue we show the region where our 
linear approach loses its validity. The excluded region is depicted 
in yellow and determined by the condition e < 2.3 X 10 22 GeV -2 . 
Finally, the density marking the beginning of structure formation 
(SF) and the dark energy (DE) density (p DE ~ 2.8 X 10~ 47 GeV 4 ) 
have also been plotted for reference. 



Model 2: f(R) = eR' 1 

We will continue our analysis with the f(R) model 
proposed in reference 27] as a dark energy candidate. 
This possibility is currently excluded, but this model is a 
simple example that help to understand the gravitational 
collapse modifications in models that provide late-time 
acceleration. 

For this model, equation (pM)) becomes: 



9'(ip) + ff(V0 csc(V0 = 



whose full solution is 



6fc 2 



sin 



cos 



13 



(34) 
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Figure 3: Analogous representation to the one shown in Fig. [T] 
which includes the GR solution, the modification given by 1 1351) and 
the sum of the two. In this figure g{ip) = 10 19 g(tf>) in order to make 
the modification observable. 
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Figure 4: Analogous to Figure [2] for model 2. The limit of the 
region depicted in yellow shows the proposed value e = — /i 4 , and 
H = 1CT 42 GeV according to 



In Figure G2 it is possible to see the behavior of the 
modification for e = — /i 4 , and /i = 10~ 42 GeV as it was 
the value originally proposed in reference (27| . 

The series expansion of (|35|) around ip = tt reads, in 
this case: 



llire(4> — it) 
8192fc 2 



(36) 



Once again the intersection between the particular so- 
lution of (fM)) and the GR solution Aq can be deter- 
mined in the ip — > ir limit, with help of equation (1311) . 
\g(ip)\ = \A G \ implies 



IItt \e\ 
2048 fc 2 ' 



(37) 



As it can be seen in Figure El the difference between 
the modified ipe and ipe, GR is not distinguishable for 
e = —/x 4 if density is higher than the standard dark 
energy density. The same result is found for e > — /i 4 



(|e| < fi 4 and negative). The situation changes for 
e < — /i 4 (|e| > /i 4 and positive). This behavior can 
be observed in Figure [H where the validity of the linear 
regime is shown to decrease for higher values of |e| and 
lower densities. As we will see in the following example, 
this is a general property of f(R) models that provide ac- 
celerated cosmologies, at least, for densities higher than 
the vacuum energy. Results in Figure |4] can be under- 
stood by estimating the correction of the collapsing time 
in the linear regime as it is determined by equation (|30l) . 
This f(R) model provides a relative difference for the 
collapsing time in GR value given by 



At 

tGR 



11 3 £ 3 7T 2 

3fc 6 2 34 



(38) 



We observe that the correction is more negligible for 
denser objects. This unexpected fact can be understood 
since GR modification to the scale factor is proportional 
to g (x e/k 2 whereas k cx po- According to this depen- 
dence, a stellar object with a higher density will suffer 
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a less important modification and vice versa. The rela- 
tive time modification is lower for denser media since the 
correction decreases with p as At/t GR oc 6 . 



Model 3: f(R) = XR 



1 + 



to be of the order of the presently observed effective cos- 
mological constant 3 . With such parameter choice, this 
model is a viable dark energy candidate. The relation 
between A and Rq in vacuum is given by Hq = XRq/6 
according to [Hf , where H is the present Hubble param- 
eter (see [29[ for recent WMAP data) and the proposed 
value for A = 0.69. 



The last f(R) model to be analyzed in the present work 
is the well-known Starobinsky model proposed in Refer- 
ence [28|. For this model, n, A > and Ro is considered 

3 It is important to remark that Ro in this context is a parameter 
of the model and not the initial scalar curvature. 



For the sake of simplicity, let us choose n = 1. In this 
case, the equation (j24|) may be rewritten as follows: 



9fcAi? o sin 3 (V0csc 4 ($) 

g'W ~ fl#) csc(V>) ^- J - 

32 (9k 2 + R 2 f 

9fcAi? 3 tan (f ) sec 4 (f ) (it* - 27k 2 sec 12 (f 
2 (9fc 2 sec 12 (|) + R 2 



{Rl + 3k 2 



= 0. 



r 



72k\R< 



o sin 4 ($J csc(-0) 
(i? 2 + 9fc 2 sec 12 ) 



Rl + 3k 2 sec 12 



(39) 



Unlike the other two cases, we are not able to find an- 
alytical solution for equation (l39l) . Thus, specific values 
for k and A parameters and Rq are required to find a nu- 
merical solution. This solution is plotted in Figure O In 
any case, equation (|39j) can be studied in the asymptotic 
limits ip — ¥ and ip — > tt. Thus, the corresponding series 
expansion of (13T)|) in the ip — > becomes 



g(iP) 9k\R (3fc 2 



R*)iP 



4(9fc 2 + ^) 2 



0. 



whose analytical solution is 



d 3k\R (R 2 - 3k 2 ) ip 2 

g(ip) = — H 5 

V; 4(9P + i? 2 ) 2 



(40) 



(41) 



Since the homogeneous equation does not depend on the 
f(R) model, the condition c\ — is also necessary in 
order to have a finite solution. When the considered 
asymptotic limit is ip — ¥ 7T, equation (|39j) approximately 
becomes 



g(iP) 9kXR (ijj - tt) 3 (3k 2 + R 2 ) 
^-tt 32(9k 2 + R 2 f 



= 0. 



(42) 



whose analytical solution is 

3kXR {ip - tt) 4 (3fc 2 



g{ip) =ci(ip-n) + 



2R 2 ) 



32 (9k 2 +R 2 Y 



(43) 



This asymptotic limit of the linear correction depends on 
a higher power of (ip — tt) than the GR solution given by 
eq. (|3"Tj) . This fact implies that we cannot estimate the 
validity of the linear regime by using the ip — ¥ tt as in the 
previous cases. The modification is more important at 
intermediates values of ip, as it can be observed in Fig- 
ure [5j The numerical results are showed in Figure 6. In 
a similar way to the second model, for A < Ao the modifi- 
cation of the collapse is always linear and not important. 
The situation is different for A > Ao , where the collapse is 
severely modified at densities closer to the vacuum one. 
We have checked numerically that denser environments 
are less affected by this gravitational model. 



V. CONCLUSIONS 

In this work we have studied the gravitational collapse 
in /(it!) gravity theories. These theories provide correc- 
tions to the field equations that modify the evolution of 
gravitational collapse with respect to the usual General 
Relativity results. In this context, viable f(R) models 
must provide similar results for the collapse times to the 
values obtained in General Relativity. In addition, col- 
lapse times must be much shorter than the age of the 
universe and long enough to allow matter cluster. 

The analyzed f(R) models present both important dif- 
ferent quantitative and qualitative behaviors when com- 
pared with General Relativity collapses. In fact, all of 
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Figure 5: The plotted lines are analogous to the ones in Fig. [T] 
and Fig. f3] but with the solution given by H39I I. In this figure 
g(ip) = 10 9 g(V>) in order to make the modification observable. 



them show a collapsing initial epoch with higher contrac- 
tion than in General Relativity. This result is expected 
since f(R) theories modify the gravitational interaction 
by the addition of a new scalar mediator. It is well-known 
that a scalar force is always attractive and can only re- 
duce the time of gravitational collapse. This result is 
interesting since observations of structures at high red- 
shift introduce some tension with the standard ACDM 
model 'Mi . and the tendency of f(R) models to increase 
the gravitational attraction at early times can alleviate 
this problem. 

Although this general behavior is shared by the three 
models analyzed throughout this investigation, they 
present significant differences when the modifications to 
the General Relativity collapse leave the linear regime. 
On the one hand, the R 2 model has a modification 
that increases with the density of the collapse object: 
(Ate/tan) oc yfp. The opposite behavior is found for the 
R^ 1 model, where this modification decreases with den- 
sity as (At c / ten) oc p~ 6 . Finally, a similar situation is 
reproduced numerically in the Starobinsky model. The 
departure from the linear collapse is able to exclude inter- 
esting parameters regions of these models that support 
late-time acceleration as seen in Figures 2, 4 and 6. 

Another relevant question is related to the physics of 
stellar objects when analyzed in the f(R) modified grav- 
ity theories frame. Although we cannot use straightfor- 
wardly the results of this analysis due to the fundamental 
role that pressure plays in the stability of these objects, 
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Figure 6: Analogous to Figure [2] for studied model 3. The limit 
of the region depicted in yellow shows the proposed value A = 0.69 
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we may get an idea of the importance of the corrections. 
Inside these objects, pressure is the same order of mag- 
nitude as density, and it is expected to introduce an im- 
portant modification into star evolution and dynamics. 
Therefore, it is enough to take into account the typi- 
cal value of the density of a neutron star, approximately 
10~ 3 GeV 4 to estimate if correction will be important. 
Although this value is 35 orders of magnitude larger than 
the dust density used above, the results do not change 
dramatically. A direct extrapolation suggests that we 
may expect even more negligible modifications for f(R) 
models that present dark energy scenarios (models 2 and 
3). In models with higher powers of the scalar curva- 
ture, the correction to General Relativity will be more 
important but still negligible (as for model 1). 
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